Theoretical prediction of perfect spin filtering at interfaces between close-packed 
surfaces of Ni or Co and graphite or graphene 
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The in-plane lattice constants of close-packed planes of fee and hep Ni and Co match that of 
graphite almost perfectly so that they share a common two dimensional reciprocal space. Their 
electronic structures are such that they overlap in this reciprocal space for one spin direction only 
allowing us to predict perfect spin filtering for interfaces between graphite and (111) fee or (0001) 
hep Ni or Co. First-principles calculations of the scattering matrix show that the spin filtering is 
quite insensitive to amounts of interface roughness and disorder which drastically influence the spin- 
filtering properties of conventional magnetic tunnel junctions or interfaces between transition metals 
and semiconductors. When a single graphene sheet is adsorbed on these open d-shell transition metal 
surfaces, its characteristic electronic structure, with topological singularities at the K points in the 
two dimensional Brillouin zone, is destroyed by the chemical bonding. Because graphene bonds only 
weakly to Cu which has no states at the Fermi energy at the K point for either spin, the electronic 
structure of graphene can be restored by dusting Ni or Co with one or a few monolayers of Cu while 
still preserving the ideal spin injection property. 

PACS numbers: 72.25.Pi,72.15.Gd,75.50.Rr 
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I. INTRODUCTION 



We recently predicted a perfect spin filtering effect 
for ultra-thin films of graphite sandwiched between 
two ferromagnetic leads. ^ This prediction emerged from 
two rapidly developing branches of condensed matter 
physics: magnetoelectronics^ and graphene electronicsi^ 
Magneto-electronics exploits the additional degree of 
freedom presented by the intrinsic spin and associated 
magnetic moment of electrons while graphene electronics 
is based upon the unique electronic properties of two- 
dimensional graphene sheets. Based on the giant magne- 
toresistance effect discovered twenty years ago,'*'^ mag- 
netoelectronics was rapidly applied to making improved 
read head sensors for hard disk recording and is a promis- 
ing technology for a new type of magnetic storage device, 
a magnetic random access memory. The giant magne- 
toresistance (GMR) effect is based on the spin depen- 
dence of the transmission through interfaces between nor- 
mal and ferromagnetic metals (FM). The effect is largest 
when the current passes through each interface in a so- 
called current-perpendicular-to-the-plane (GPP) measur- 
ing configuration but the absolute resistance of metallic 
junctions is too small for practical applications and the 
current-in-plane (GIP) configuration with a much smaller 
MR is what is used in practice. Replacing the non- 
magnetic metal spacer with a semiconductor^ or insulator 
(I), such as Al2 03'''^ results in spin-dependent tunneling 
and much larger resistances are obtained with FM|I|FM 
magnetic tunnel junctions (MTJs). Substantial progress 



has been made in increasing the tunneling MR effect 
by replacing the amorphous AI2O3 insulator with crys- 
talline MgO,^'--'^^ Though there is a relatively large lattice 
mismatch of 3.8% between Fe and MgO, the tunneling 
magnetoresistance (TMR) in Fe|MgO|Fe junctions has 
been reported to reach values as high as 180% at room 
temperature.— Low temperature values as high as 1010% 
have been reported for FeCoB|MgO|FeGoB MTJsJ^ii^ 
The sensitivity of TMR (and spin injection) to details 
of interface structur e — .i -'^ makes it difficult to close the 
quantitative gap between theory and experiment so it is 
important for our understanding of TMR to be able to 
prepare interfaces where disorder does not dominate the 
spin filtering properties. This remains a challenge due 
to the high reactivity of the open-shell transition metal 
(TM) ferromagnets Fe, Go, and Ni with typical semicon- 
ductors and insulators. 

With this in mind, we wish to draw attention to a quite 
different material system in which a thin graphite film is 
sandwiched between two ferromagnetic leads. Graphite 
is the ground state of carbon and as one of the most 
important elemental materials, its electronic structure 
has been studied in considerable detail. It consists 
of weakly interacting sheets of carbon atoms strongly 
bonded in a very characteristic honeycomb structure. Be- 
cause of the weak interaction between these "graphene" 
or "monolayer graphite" sheets, the electronic structure 
of graphite is usually discussed in two steps: first, in 
terms of the electronic structure of a single sp^ -bonded 
sheet, followed by consideration of the interaction be- 
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FIG. 1: Fermi surface projections onto close-packed planes for: (a) fee Cu; (c) majority- and (d) minority-spin fee Ni (111); 
(e) majority- and (f) minority-spin fee Co (111); (g) majority- and (h) minority-spin hep Ni (0001); (i) majority- and (j) 
minority-spin hep Co (0001). For graphene and graphite, surfaces of constant energy are centred on the K point of the two 
dimensional interface Brillouin zone (b). The number of Fermi surface sheets is given by the colour bar. 



tween sheets.— li^i^ From these early, and many subse- 
quent studies, it is known that graphene is a "zero-gap 
semiconductor" or a semimetal in which the Fermi sur- 
face is a point at the "K" point in the two-dimensional 
reciprocal space. The physical properties associated 
with this peculiar electronic structure have been stud- 
ied theoretically in considerable detail, in particular in 
the context of carbon nanotubes which can be consid- 
ered as rolled-up graphene sheets jiS, With the very re- 
cent discovery and development of an exceptionally sim- 
ple procedure for preparing single and multiple graphene 
sheets, micromechanical cleavage,— it has became pos- 
sible to probe these predictions experimentally. Sin- 
gle sheets of graphene turn out to have a very high 
mobility^ ^ that manifests itself in a variety of spectac- 
ular transport phenomena such as a minimum conduc- 
tivity, anomalous quantum Hall effect (QHE),^^'^'^ bipo- 
lar supercurrent24 and room-temperature QHEJii Spin 
injection into graphene using ferromagnetic electrodes 



TABLE I: Lattice constants of Co, Ni, Cu, and graphite, 
ahox = atcc/\/2. Equilibrium separation do for a single 
graphene sheet on top of the graphite (0001) and Co, Ni or 
Cu fee (111) surfaces as calculated within the framework of 
the DFT-LDA using the in-plane lattice constant Ohcx ~ 2.46 

A. 





Graphite 


Co 


Ni 


Cu 


a-f (A) 

<r (A) 

alg^ (A) 
do (A) 


2.46 
2.45 
3.32 


3.544" 
2.506 
2.42 
2.04 


3.524" 
2.492 
2.42 
2.03 


3.615" 
2.556 
2.49 
3.18 





has already been realized . ^^'^^ The weak spin-orbit in- 
teraction implied by the low atomic number of carbon 
should translate into very long intrinsic spin-flip scatter- 
ing lengths, a very desirable property in the field of spin 
electronics or "spintronics" , which aims to combine tradi- 
tional semiconductor-based electronics with control over 
spin degrees of freedom. However, the room tempera- 
ture two-terminal MR effect of ^ 10% observed in lat- 
eral, current-in-planc (CIP) graphene-based devices with 
soft permalloy leads is still rather small^^S 

Instead of a CIP geometry, we consider a CPP 
TM|Gr|TM (111) junction, where TM is a close-packed 
surface of fee or hep Ni or Co and Gr is graphite (or n 
sheets of graphene, Gr„). We argue that such a junction 
should work as a perfect spin filter. The essence of the 
argument is given by Table U and Fig. [1] According to 
Table|Tl the surface lattice constants of (111) Ni, Co and 
Cu match the in-plane lattice constants of graphene and 
graphite almost perfectly. The lattice mismatch of 1.3% 
at the Ni(lll)|Gr interface is, in fact, one of the small- 
est for the magnetic junctions that have been studied so 
far. This small lattice mismatch suggests that epitaxial 
TM|Gr|TM junctions might be realized experimentally, 
for example using chemical vapor depositioni^i^SiSi As- 
suming perfect lattice matching at the TM|Gr interface, 
it is possible to directly compare the Fermi surface pro- 
jection of graphite with the projections of the Fermi sur- 
faces (FS) of fee Cu and of fee and hep Ni and Co onto 
close-packed planes, see Fig. [T] 

The Fermi surface of graphene is a point at the high- 
symmetry K point in reciprocal space. The Fermi sur- 
faces of graphite and of doped graphene are centred on 
this point and close to it. Figure[T]shows that there are no 
majority spin states for Ni and Co close to the K point 
whereas minority spin states exist (almost) everywhere 
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in the surface BZ. Only the minority spin channel should 
then contribute to transmission from a close-packed TM 
surface into graphite. In a TM|Gr|TM junction, electrons 
in other regions of reciprocal space on the left electrode 
would have to tunnel through graphite to reach the right 
electrode. If the graphite film is taken thick enough to 
suppress tunneling, majority spin conductance will be 
quenched and only minority spin conductance through 
the graphite will survive i.e. perfect spin filtering will 
occur when the magnetizations are aligned in parallel 
(P). For antiparallel (AP) alignment, the conductance 
will vanish. 

In this paper, we wish to study the effectiveness of 
this spin filtering quantitatively: how it depends on the 
thickness of the graphite film, the geometry of the clean 
metal-graphite interface, interface roughness and disor- 
der, and lattice mismatch. While we will be mainly con- 
cerned with the CPP geometry, we will also comment 
on the applicability of some of our conclusions to the 
CIP geometry. The paper is organized as follows. In 
Sec. In] we give a brief description of the computational 
method and outline the most important technical de- 
tails of the calculations. The transport formalism we 
use is based upon a very efficient minimal basis of tight- 
binding mufhn tin orbitals (TB-MTO) in combination 
with the atomic spheres approximation (ASA).^^ While 
the ASA works well for close-packed structures, some 
care is needed in using it for very open structures like 
that of graphite. In Sec. IIIII we therefore benchmark the 
electronic structures calculated using the TB-MTO- ASA 
method with those obtained from plane-wave pseudopo- 
tential calculations. Section IIVI contains the results of 
spin-dependent electron transport calculations for spec- 
ular interfaces (ideal junction) as well as for junctions 
with interface roughness and alloy disorder. A summary 
is given and some conclusions drawn in Sec. |Vl 



II. COMPUTATIONAL METHOD 

The starting point for our study is an atomic struc- 
ture calculated by minimizing the total energy within 
the local spin density approximation (LSDA) of density 
functional theory (DFT). This was done using a plane- 
wave pseudopotential (PWP) method based upon pro- 
jector augmented wave (PAW) pseudopotentials'^'^ as im- 
plemented in the VASP program i^^i^^i^ The interaction 
between graphite and the TM surface is modelled us- 
ing a repeated slab geometry of six metal layers with a 
graphene sheet on top and a vacuum thickness of 12 
A. To avoid interactions between periodic images of the 
slab, a dipole correction is applied.^'' The surface Bril- 
louin zone (SBZ) was sampled with a 36 x 36 k-point grid 
and the SBZ integrals carried out with the tetrahedron 
integration scheme.— A plane wave kinetic energy cutoff 
of 400 eV was used. The plane-wave pseudopotential cal- 
culations yield energy band structures, charge transfers, 
binding energies and work functions for single TM|Gr 



interfacesi^ The equilibrium distances do between the 
graphene sheet and the TM surfaces are summarized in 
Table H 



The equilibrium geometries are used as input for self- 
consistent TB linearized MTO (TB-LMTO)^^ calcula- 
tions for the TM|Gr„|TM junction. The resulting Kohn- 
Sham potentials are used to calculate spin-dependent 
transmission probabilities through the TM|Gr„|TM junc- 
tion using a TB-MTO wave-function matching^i^ 
scheme42ii^i^ To do this, the junction is divided into 
three parts consisting of a scattering region sandwiched 
between semi-infinite left and right leads, all of which are 
divided into layers that are periodic in the lateral direc- 
tion. The leads are assumed to be ideal periodic crystals 
in which the electron states (modes) are wave functions 
with Bloch translational symmetry. By making use of 
its Bloch symmetry, a semi-infinite lead can be repre- 
sented as an energy-dependent non-Hermitian potential 
on the boundary of the scattering region so that the in- 
finite system is made finite. According to the Landauer- 
Biittiker formalism of transport, the conductance can be 
calculated by summing up all the probabilities for trans- 
mitting an electron from the electron modes in the left 
lead through the junction into electron modes in the right 



The effect of various types of disorder on the trans- 
mission can be studied using the same formalism and 
computer codes by modelling the disorder within large 
lateral supercells^i^ and averaging over many configu- 
rations of disorder generated by choosing positions of im- 
purity atoms or imperfections randomly. We study three 
types of disorder: interface roughness, interface alloying 
and lattice mismatch. In the first two cases, averaging is 
performed over a minimum of ten configurations of dis- 
order. To model interface roughness, some surface atoms 
are removed (replaced by "empty spheres" with nuclear 
charges that are zero in the ASA) and the ASA potentials 
are calculated self-consistently using a layer version**^ of 
the coherent potential approximation (CPA) 4^ The ef- 
fect of interface alloying which might occur if deposition 
of a thin layer of Cu on Ni or Co ("dusting") leads to 
intermixing is modelled in a similar fashion. Thirdly, the 
small lattice mismatch between graphite and TM is mod- 
elled by "cutting and pasting" AS potentials from self- 
consistent calculations for TM|Gr„|TM junctions with 
two different in-plane lattice constants. The two systems 
are then combined using a supercell whose size is deter- 
mined by the lattice mismatch. For self-consistent TB- 
LMTO-ASA calculations, the BZ of lateral supercells is 
sampled with a density roughly corresponding to a 24 x 24 
k-point grid for a 1 x 1 interface unit cell. To converge the 
conductance, denser grids containing 800 x 800, 20 x 20 
and 8x8 k-points are used for 1x1 (ideal junction), 
5x5 and 20 x 20 lateral supercells, respectively. 
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FIG. 2: Top and side perspective views (top and bottom 
panels) of graphite where the potential is represented in the 
atomic spheres approximation using additional, empty atomic 
spheres. Model I (left) contains 32 empty spheres in a unit 
cell containing 4 carbon atoms (red spheres). Model II (right), 
contains just 4 empty spheres. For model I, gray, green, blue 
and yellow spheres display the positions of the empty spheres 
El, E2, E3 and E4, respectively. For model II, there is just 
one type of empty sphere (green). 



III. GEOMETRY AND ELECTRONIC 
STRUCTURE OF TMjGrnlTM 

In this section we describe in more detail how the elec- 
tronic structure of TM|Gr„|TM junctions for TM=Cu, 
Ni or Co is calculated. These close-packed metals can 
be grown with ABC stacking in the (111) direction (fee), 
or with AB stacking in the (0001) direction (hep). We 
neglect the small lattice mismatch of 1.3%, 1.9% and 
3.9% for the Ni|Gr, Co|Gr, and Cu|Gr interfaces, respec- 
tively, and assume the junction in-plane lattice constant 
to be equal to that of graphite, acr = 2.46 A. In the 
atomic spheres approximation, the atomic sphere radii 
of Ni, Co and Cu are then ttm = 2.574 a.u. The ASA 
works well for transition metals like Co, Ni or Cu which 
have close-packed structures. For materials like graphite 
which has a very open structure with an in-plane lat- 
tice constant acr = 2.46 A, and an out-of-plane lat- 
tice constant cqi- = 6.7 A, the unmodified ASA is not 
sufficient. Fortunately, a reasonable description of the 
crystal potential can be obtained by packing the intersti- 
tial space with empty spheres4^ This procedure should 
satisfy the following criteria: (i) the total volume of all 
atomic spheres has to be equal to the volume of the en- 
tire system (space filling), and the (ii) overlap between 
the atomic spheres should be as small as possible. 



TABLE II: Wyckoff symbols, standardized position parame- 
ters and atomic sphere radii for carbon atoms, C, and empty 
spheres, E (with nuclear charge Z—Q), for two structural mod- 
els of graphite with space group 7?|^ (PGa/mmc) Reflsol. 
Model I contains four different types of empty sphere: Ei, 
E2, E3, E4; model II only one, E. 



Model 


Atom 


Wyckoff 


position 


radius 






position 


parameters 


fa u ) 


I 


Ci 


2b 




1.56 




C2 


2c 




1.56 




El 


2a 




1.4 




E2 


2d 




1.6 




E3 


4f 


z=0.5 


1.4 




E4 


241 


x=l/3, y=0, z=0.38 


0.9 


II 


Ci 


2b 




1.56 




C2 


2c 




1.56 




E 


4f 


z=0.4 


2.18 



A. Graphite and graphene 

To see how this procedure works in practice, we bench- 
mark the TB-MTO-ASA band structure of graphite 
against the "exact" band structure calculated with the 
PWP method. To preserve the graphite IJg^ (P63/mmc) 
space group symmetry,— the positions of the atomic 
spheres are chosen at Wyckoff positions. There are twelve 
different Wyckoff positions consistent with Df/^ symme- 
try and the best choice of empty spheres is not immedi- 
ately obvious. We construct two models that describe the 
band structure close to the Fermi energy well compared 
to the PWP results; this is what is most relevant for 
studying transport in the linear response regime. Model 
I with 32 empty spheres per unit cell and model II with 
only 4 empty spheres per unit cell both preserve the sym- 
metry of graphite within the ASA. The crystal structures 
of graphite packed with empty spheres according to these 
two models is shown schematically in Fig. [51 Note that 
not all the empty spheres in a unit cell are shown in the 
figure. The Wyckoff labels, atomic sphere coordinates 
and radii are given in Table [TTl Figure [3] shows the band 




FIG. 3: (Color online) Band structure of graphite for model 
I (on the left), and model II (on the right). Gray (green) 
dots and black lines correspond to band structures calculated 
using the PWP and TB-MTO-ASA methods, respectively. 
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structure of graphite obtained with the TB-MTO-ASA 
for models I and II compared to the "exact" PWP band 
structure. Both models are seen to describe the graphite 
TT bands around the Fermi energy very well. Model I pro- 
vides a very good description of the bands within ±2 eV 
of the Fermi energy, while the smaller basis model II is 
quite good within ±1 eV. At the cost of including many 
more empty spheres, model I provides a better descrip- 
tion of the crystal potential between the graphene planes 
than model II. For this reason we use model I to study 
the transport properties of ideal junctions, junctions with 
interface roughness and alloy disorder. To be able to han- 
dle the large 20 x 20 lateral supercells needed to model 
a lattice mismatch of 5% at the TM|Gr interface, we use 
model II. 



B. Graphene on Ni(lll) substrate 

The next step is to put a monolayer of graphite 
(graphene) on top of the Ni(lll) substrate at a distance 
do from the metal surface. From our studies of the ener- 
getics of graphene on TM(lll), we found^i^ that the low- 
est energy configuration (with 3m symmetry) for TM=Ni 
or Co corresponds to an "AC" configuration in which one 
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FIG. 4; (Color online) "AC" model of TMjGr„|TM structure 
for (a) odd and (b) even numbers of graphene sheets. Car- 
bon atoms are represented by small dark (red) spheres, TM 
atoms by larger gray spheres. The configuration shown in 
(a) is a ciCi configuration with the carbon atom labelled ci 
above an "A" site surface layer TM atom of the top and the 
bottom electrodes. The other carbon atom, C2, is above a 
third layer TM atom on a "C" site. An equivalent C2C2 con- 
figuration in which the C2 atoms are on top of "A" site TM 
atoms can be realized by rotating the top and bottom elec- 
trodes by 180° about a vertical axis through the second layer 
"B" sites; this effectively interchanges ci and C2. Two other 
equivalent configurations C1C2 and C2C1 can be realized in an 
analogous fashion by rotating either the top or the bottom 
electrode through 180°. For two sheets of graphene stacked 
as in graphite, a C2C2 configuration is sketched in (b). In- 
terlayer distance is indicated as do and c/2 is the distance 
separating two neighbouring graphene sheets. 



carbon atom is positioned on top of a surface TM atom 
(an "A" site) while the second carbon atom is situated 
above a third layer TM atom (a "C" site), where A and 
C refer to the ABC stacking of fee close-packed planes, 
see Fig. |4l This is in agreement with another recent first- 
principles calculations^ as well as with experiment a^°i'^^ 
for graphene on the Ni(lll) surface. The electronic struc- 
ture of a single graphene sheet will depend on do and the 
details of such graphene-metallic substrate contacts can 
be expected to play an important role in current-in-plane 
(CIP) devicesj^iSl For the less strongly bound BC config- 
uration of Gr on Ni, the equilibrium separation is rather 
large, do ~ 3.3 A and the characteristic band structure 
of an isolated graphene sheet is clearly recognizable; see 
Fig. O For the lowest energy AC configuration, the in- 
teraction between the graphene sheet and Ni surface is 
much stronger, a gap is opened in the graphene derived 
Pz bands and at the Fermi energy there are no graphene 
states at the K-point in reciprocal space for the minority 
spin channel. This may complicate efficient spin injection 
into graphene in lateral, CIP devices^SS, 

The band structure calculated with the TB-MTO-AS 
approximation for the AC configuration is shown in the 
bottom panel of Fig.[5]and is seen to describe graphene on 
Ni(lll) qualitatively quite well. However, the splitting of 
the graphene bands, which arises because the two carbons 
atoms are no longer equivalent when one is above a top 
layer A site Ni atom and the other is above a third layer 
C site Ni atom, is somewhat larger than that resulting 
from the PWP calculation. 



C. Ni|Gr„|Ni(lll) junction 

The transmission of electrons through a TM|Gr„|TM 
junction will obviously depend on the geometry of the 
metal-graphite contacts. Rather than carrying out a 
total energy minimization explicitly for every different 
value of n, we assume that the weak interaction between 
graphene sheets will not infiuence the stronger TM|Gr 
interaction and construct the junction using the "AC" 
configuration and the equilibrium separation do = 2.03 
A for each interface, as shown in Fig. 21 The interstitial 
space at the TM | Gr interfaces is filled with empty spheres 
using a procedure analogous to that described for bulk 
graphite 1^ 

Because the two carbon atoms Ci and C2 in the 
graphene unit cell are equivalent, either of them can be 
positioned above a surface Ni atom on an A site with the 
other on the C site in an "AC" configuration, without 
changing the total energy. Since this can be done for 
each TM|Gr interface separately, four different configu- 
rations of the TM|Gr|TM junction can be constructed 
by rotating one or both electrodes through 180° about a 
vertical axis through the second layer B sites which inter- 
changes electrode A and C sites in Fig. [H We label these 
four different configurations CiCi, C1C2, C2C1 and C2C2 in 
terms of the carbon atoms which are bonded to A site 
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FIG. 6: Energy band structures of an ideal NielGr (111) mul- 
tilayer with 6 layers of fee Ni sandwiching a single graphene 
sheet in a ciC2 configuration, for majority (left panels) and 
minority (right panels) spin channels. Plane wave pseudopo- 
tential calculations are shown on top (dotted lines), the TB- 
LMTO-AS A results on the bottom (soUd Unes) . 



FIG. 5: (Color online) The results of PWP (top and mid- 
dle rows) and TB-MTO-ASA (bottom row) calculations of 
majority (left panels) and minority (right panels) spin band 
structures (green) of single graphene layers absorbed on both 
sides of a 13 layer (111) Ni slab for a BC configuration with 
do = 3.3 A, (top) and an AC configuration with do = 2.0 A 
(middle and bottom). The bands are replotted and superim- 
posed in black using the carbon pz character as a weighting 
factor. The Fermi energy is indicated by the horizontal dashed 
line. 



TM atoms. For more than one graphene sheet, the sec- 
ond sheet breaks the symmetry between the Ci and C2 
atoms. While we have not checked this explicitly, we ex- 
pect the corresponding energy difference to be small and 
neglect it. 

In Figure [5] we saw that the graphene tt states inter- 
acted strongly with the nickel surface in the minimum en- 
ergy "AC" configuration. The interaction with the metal 
substrate made the Ci and C2 carbon atoms inequivalent 
and led to the opening of an energy gap in the graphene 
TT bands. Having constructed an interface geometry, we 
study the band structure of the Ni| graphene |Ni (111) 
junction as a function of k|| , the two dimensional Bloch 
vector, modelling it as a Nisj graphene jNis junction re- 
peated periodically in the (111) direction (which is equiv- 
alent to a NigjGri multilayer). The bands in the top pan- 
els of Fig. [6] were calculated using the benchmark plane- 
wave pseudopotential (PWP) method, those in the bot- 
tom panels with the TB-LMTO-ASA. We see that the 
Ni-related bands are described well by the TB-LMTO- 
ASA - as might be expected since the ASA is known to 
work well for close-packed solids. The second thing we see 
is that there is no gap in the graphene tt bands. This is 



because the C1C2 configuration used in the calculation has 
inversion symmetry and the equivalence between the two 
carbon atoms is restored; see Fig. Ufa). The third point 
to be made is that the charge transfer from graphene 
(work function: 4.5 eV) to Ni (work function: 5.5 eV) 
and strong chemisorption leads to the formation of a po- 
tential step at the interface and a significant shift of the 
graphene tt bands with respect to the Fermi level^ which 
is pinned at that of bulk Ni. We find similar results for 
Co I Gr I Co ( 1 1 1 ) and Co I Gr I Co(000 1 ) junctions . There is a 
difference between the position of the graphene 7r-derived 
bands, most noticeably at the K point, in the PWP and 
TB-LMTO-ASA band structures shown in Fig.[6]for both 
spin channels. It appears that the interface dipole is not 
accurately described by the ASA. From the point of view 
of describing transmission of electrons through this junc- 
tion, the electronic band structure is the most important 
measure of the quality of our basis, description of the po- 
tential, etc., so this discrepancy will most certainly have 
quantitative implications. Fortunately, our most impor- 
tant conclusions will be qualitative and will not depend 
on this aspect of the electronic structure. 



IV. ELECTRON TRANSPORT THROUGH A 
FM|GrnlFM JUNCTION 

Using the geometries and potentials described above, 
we proceed to study the spin-dependent transmission 
through ideal Ni|Gr|Ni junctions in the CPP geometry as 
a function of the thickness of the graphite spacer layer. 
We then discuss how interface roughness, alloy disorder 
and the lattice mismatch between graphite and the sub- 
strate affect the spin-filtering properties of the junctions 
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using large lateral supercells to model the various types 
of disorder. Because very similar results are found for all 
the TMs shown in Fig. [l] we focus on foe Ni as a sub- 
strate because it has the smallest lattice mismatch with 
graphite and graphene has been successfully grown on Ni 
using chemical vapour deposition i^^'^°i'^^ 



A. Specular interface 

The spin-dependent transmission through Ni|Gr„|Ni 
(111) junctions is shown in Fig.[7]for parallel (P) and an- 
tiparallel (AP) orientations of the magnetization in the 
nickel leads, in the form of the conductances Gp and 
G^p with cr = min, maj. All the conductance values 
are averaged over the four interface configurations of the 
Ni|Gr„|Ni junction which are consistent with AC con- 
figurations of the Ni|Gr (111) interface. Gp^'' and GJp 
are strongly attenuated, while Gp™ saturates to an n- 
independent value. The magnetoresistance (MR) defined 
as 

MR = ' X 100% ^ ^''^ ^^"^ X 100%. (1) 
Rap Gp 

rapidly approaches its maximum possible value of 100%, 
as shown in the right inset in Fig. [71 This pessimistic 
definition of MR is more convenient here because Gap 
vanishes for large n. It is usually the optimistic version, 
that approaches 10^^ % in our calculations but does not 
saturate, that is quoted. ■^'i^'^'^ The left inset in Fig. [7] 
shows how the conductance depends on the particular 
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FIG. 7: Conductances Gp'" (V), Gp^-^ (a), and Gap (x) av- 
eraged over the four configurations ciCi, ciC2, C2C1 and C2C2 of 
a Ni|Gr„|Ni junction as a function of the number of graphene 
monolayers n for ideal junctions. Right inset: magnetoresis- 
tance MR as a function of n. Left inset: minority parallel 
conductance Gp'" (V) given for four different configurations. 
The points which are circled and connected with a dashed line 
are the values which were shown in Ref. [J. 



configuration of the junction. The minority spin conduc- 
tance in the parallel configuration, which dominates the 
magnetoresistance behaviour, is highest for the ciCi con- 
figuration with an asymptotic value of Gp"" ~ 10^^ Gq. 
This is approximately an order of magnitude larger than 
Gp"" for the C2C1, C1C2 and C2C2 configurations. The C1C2 
and C2C1 configurations are equivalent so the correspond- 
ing values of Gp"" should be identical. The small dif- 
ferences between these two configurations which can be 
seen in the figure are an indication of the overall accuracy 
of the numerical calculation. The points which are cir- 
cled and connected with a dashed line are the oscillating 
values which were shown in Ref. [l]. 

To demonstrate that spin-filtering occurs due to high 
transmission of minority spin electrons around the K 
point, we plot the majority- and minority-spin transmis- 
sion for the P configuration as a function of k|| for two 
graphite films of different thickness in Figure [S] A single 
sheet of graphene (a monolayer of graphite) is essentially 
transparent with a conductance of order Go in both spin 
channels. In the minority spin channel, the transmission 
is very low or vanishes close to F and close to K along 
the high symmetry F-K line, in spite of there being one 
or more sheets of Fermi surface in these regions of recip- 
rocal space. This is a clear indication of the importance 
of matrix element effects: selection rules resulting from 
the incompatibility of wave functions on either side of the 
interfacCfi? 

The majority transmission must be zero around F and 
around the K point because there are no states there in 
the Ni leads. For thicker graphite, the only contribution 
to the majority-spin conductance comes from tunneling 
through graphite in regions of the 2D-BZ where there 



1 ML 5 ML 




FIG. 8: Transmission as a function of the transverse crystal 
momentum ky in the two dimensional interface BZ for a C1C2 
configuration of an ideal Ni|Gr„ |Ni (111) junction in a parallel 
state, (a) and (b) are for a single graphene sheet, n = 1; (c) 
and (d) are for n — 5; (e) shows the minority spin transmission 
in a small circle of radius r = 0.057 (27r/aGr) around the K 
point for 5 ML of graphite on an enlarged scale. 
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are Ni states and the gap between graphite bonding and 
antibonding tt states is small. This occurs close to the 
M point 1^ see Fig. |3l Because the gap decreases going 
from M to K, the transmission increases in this direction. 
At the edge of the Fermi surface projection, the velocity 
of the Bloch electrons in the leads is zero so that the 
maximum transmission occurs just on the M side of these 
edges. 

The total minority transmission consist of two con- 
tributions. On the one hand there is a tunneling con- 
tribution from throughout the 2D-BZ which, depending 
on the particular k|| point, is determined by the gap in 
graphite as well as by the compatibility of the symme- 
tries, at that point, of the wave functions in Ni and in 
graphite. On the other hand there is a large transmis- 
sion from the neighbourhood of the K point coming from 
the Bloch states there in graphite. Once these have cou- 
pled to available states in Ni, this contribution does not 
change much as more layers of graphite are added. Per- 
fect spin-filtering (100% magnetoresistance) occurs when 
the tunneling contributions are essentially quenched com- 
pared to the minority spin K point contribution. For four 
MLs of graphite the polarization is within a percent of 
100% and for five MLs it is for all intents and purposes 
complete. The only discernible transmission in Figs, ^c- 
e) is found close to the K point. Magnification of this 
region in Fig. [HJe) shows a certain amount of structure 
in the transmission. This can be explained in terms of 
the multiple sheets of Ni minority spin Fermi surface in 
the vicinity of K (Fig. [T]) and the small but finite dis- 
persion of the graphite bands perpendicular to the basal 
planci^ The transmission is seen to have the threefold 
symmetry of the junction. 

The spin-filtering does not depend on details of how 
graphite is bonded to the ferromagnetic leads as long as 
the translational symmetry parallel to the metal-graphite 
interfaces is preserved. We have verified this by perform- 
ing explicit calculations (results not shown here) for junc- 
tions in the "AB" and "BC" configurations with different 
metal-graphite separations d. 



B. 



Ni|Cu,„|Gr„|Cu^|Ni (111) 



In Section IIIIl we saw that the electronic structure 
of a sheet of graphene depends strongly on its separa- 
tion from the underlying TM substrate. For Co and Ni, 
equilibrium separations of the order of 2.0 A were cal- 
culated for the lowest energy AC configuration (see Ta- 
ble m , the interaction was strong and the characteristic 
linear dispersion of the graphene electronic structure was 
destroyed, Fig. [S] For a separation of 3.3 A, the small 
residual interaction does not destroy the linear disper- 
sion. Unlike Co and Ni, Cu interacts only weakly with 
graphene, there is only a small energy difference between 
the "asymmetric" AC configuration with do = 3.3 A and 
the slightly more weakly bound "symmetric" BC config- 
uration with do = 3.4 A, and bonding to Cu preserves 



the characteristic graphene electronic structure, opening 
up only a very small gap of about 10 meV at the Dirac 
point 

Should it be desirable to avoid forming a strong bond 
between graphite and the TM electrode, then it should 
be a simple matter of depositing one or a few layers of 
Cu on e.g. Ni. Such a thin layer of Cu will adopt the 
in-plane lattice constant of Ni and graphite will bind to 
it weakly so that the electronic structure of the first layer 
of graphite will be only weakly perturbed. Because Cu 
oxidizes less readily than Ni or Co, it may be used as a 
protective layer. Cu has no states at or around the K 
point for either spin channel (Fig. [T]) so it will simply 
attenuate the conductance of the minority spin channel 
at the K point. This is demonstrated in Fig [5] where 
the magnetoresistance of a Ni|Cu„i|Gr„|Cu„i|Ni junction 
is shown as a function of the number m of layers of Cu 
when there are 5 MLs and 7 MLs of graphite. As the 
thickness of Cu is increased reducing the transmission 
of the minority-spin K point channel, the MR decreases. 
The reduction of the MR can be compensated by increas- 
ing the thickness of graphite. These conclusions are con- 
sistent with the qualitative conclusions drawn above in 
connection with Fig[T](a). 

Although the linear dispersion of the graphene bands is 
essentially unchanged by adsorption on Cu, application of 
an in-plane bias will destroy the translational symmetry 
parallel to the interface upon which our considerations 
have been based. The finite lateral size of a Ni|Cu elec- 
trode will also break the translational symmetry in a CIP 
measuring configuration and edge effects may destroy the 
spin-injection properties. 
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FIG. 9: Magnetoresistance as a function of the number of Cu 
monolayers on both left and right Ni leads in case of 5 ML 
(dashed line) and 7 ML (solid line) of graphene. 
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C. Effect of disorder 

100 

1. Lattice Mismatch 



So far, we have assumed TM and graphite lattices 
which are commensurate in-plane. In practice there is 
a lattice mismatch with graphite of 1.3% for Ni, 1.9% 
for Co and 3.9% for Cu which immediately poses the 
question of how this will affect the perfect spin-filtering. 
While lattice mismatch between lattices with lattice con- 
stants oi and 02 can in principle be treated by us- 
ing ni units of lattice 1 and n2 units of lattice 2 with 
nifli = n2a2, in practice we cannot perform calculations 
for systems with n much larger than 20 which limits us to 
treating a large lattice mismatch of 5%. To put an upper 
limit on the effect of a 1.3 — 1.9% lattice mismatch, we 
performed calculations for a NijGrs |Ni junction matching 
19 X 19 unit cells of Ni in-plane to 20 x 20 unit cells of 
graphite. The effect of this 5% lattice mismatch was to 
reduce the (pessimistic) magnetoresistance from 100% to 
90% (or ~ 900% in the optimistic definition). We con- 
clude that the actual Ni|Gr mismatch of 1.3% should not 
be a serious limiting factor in practice. 

2. Interface Roughness 

Incommensurability is not the only factor that might 
reduce the magnetoresistance. Preparing atomically per- 
fect interfaces is not possible and raises the question of 
how sensitive the perfect spin-filtering will be to inter- 
face roughness or disorder. Our studies of spin injection 
in Ref. [11 and TMR in Ref. suggest they may be 
very important and can even dominate the spin trans- 
port properties. 

The simplest way to prepare a CPP Ni|Gr|Ni junction 
would presumably be to begin with a (111) oriented Ni 
or Co crystal characterized on an atomic scale by STM or 
AFM, grow the required number of layers of graphene by 
e.g. chemical vapour deposition22ii2£i^ and after char- 
acterization of the graphene layers to then deposit the 
second Ni electrode. To prepare a CIP junction, we en- 
visage a procedure in which thin graphite layers are pre- 
pared by micromechanical cleavage of bulk graphite onto 
a Si02 covered Si wafer— into which TM (Ni or Co) elec- 
trodes have been embedded. We assume that the (111) 
electrodes can be prepared in ultrahigh vacuum and char- 
acterized on an atomic scale and that the surfaces are flat 
and defect free. Layers of graphene are peeled away until 
the desired value of n is reached. 

Assuming it will be possible to realize one essentially 
perfect interface, we have studied the effect of roughness 
at the second interface, assuming it is prepared by evapo- 
ration or some similar method. The graphite is assumed 
to be atomically perfect and all of the roughness occurs 
in the metal interface layer. We model this roughness 
as in Ref. [l3 by removing a certain percentage of the 
top layer atoms. The atomic sphere potentials are cal- 
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number of graphene monolayers 

FIG. 10: (Color online) Magnetoresistance as a function 
of n for: ideal junctions (circles); Ni|Grn|Cu5oNi5oiNi junc- 
tions where the surface layer is a disordered alloy (diamonds) ; 
Ni|Gr„|Ni junctions where the top layer of one of the elec- 
trodes is rough with only half of the top layer sites occupied 
(squares). For the rough surface layer, the error bars indicate 
the spread of MR obtained for different configurations. In- 
set: schematic representation of Ni|Gr„|Ni junction with alloy 
disorder (roughness) at the right Ni|Gr interface. Ni atoms 
are given by large gray spheres while Cu (missing) atoms in 
the case of alloy disorder (roughness) are given by large dark 
(blue) spheres. Positions of carbon atoms are represented by 
small dark (red) spheres. 



culated using the layer versionsi of the coherent poten- 
tial approximation (CPA).^^ The CPA AS potentials are 
then distributed at random with the appropriate concen- 
tration in 5 X 5 lateral supercells and the transmission 
is calculated in a CPP geometry for a number of such 
randomly generated configurations. The effect on the 
magnetoresistance of removing half a monolayer of Ni is 
shown in Fig. \W\ as a function of the number of graphite 
layers. 50% roughness at one interface is seen to reduce 
the 100% magnetoresistance to about 70% 230% op- 
timistic). 



3. Interface Disorder 

The last type of disorder we consider is a layer of in- 
terface alloy. We imagine that depositing a layer of Cu 
on Ni to prevent graphite bonding to the Ni has led to a 
layer of Ni and Cu mixing. In a worst case scenario, we 
assume all of the disorder is in the surface layer and as- 
sume this to be a MsqCuso random alloy. The potentials 
are once again calculated self-consistently using the layer 
CPA and the transmission calculated as for roughness. 
The effect on a monolayer of CuNi alloy is to reduce the 
MR to 90% (900% in the optimistic definition) for a thick 
graphite film, as shown in Fig. 1101 These results indicate 



that the momentum transfer induced by the scattering 
due to imperfections is insufficient to bridge the large gap 
about the K point in the majority spin FS projections. 

Ideahy, we should avoid interface roughness and dis- 
order altogether. Since metal surfaces can be prepared 
with very little disorder, what is required is to be able 
to perform micromechanical cleavage on a metal surface 
rather than on Si02. If this were possible, two essentially 
perfect TM|Gr interfaces could perhaps be joined using 
a method analogous to vacuum bonding^^ Alternatively, 
since graphite has a large c-axis resistivity^^ it may only 
be necessary to prepare one near-perfect Ni| graphite in- 
terface. If the graphite layer is sufficiently thick, then 
it should be possible to achieve 100% spin accumulation 
in a high resistivity material making it suitable for in- 
jecting spins into semiconductors.^*^ Because carbon is so 
light, spin-flip scattering arising from spin-orbit interac- 
tion should be negligible. 
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L R 

FIG. 11: Schematic figure of a TM|Gr|TM CIP junction in 
which the electric field is forced to be essentially perpendicular 
to the TM|graphene interface. The dashed shaded box indi- 
cates an alternative configuration with the right-hand elec- 
trode on top of the graphene. 



V. DISCUSSION AND CONCLUSIONS 

Motivated by the recent progress in preparing and 
manipulating discrete, essentially atomically perfect 
graphene layers, we have used parameter-free, materials- 
specific electronic structure calculations to explore the 
spin transport properties of a novel TM| graphite system. 
Perfect spin-filtering is predicted for ideal TM|Gr„|TM 
junctions with TM = Co or Ni in both fee and hep crystal 
structures. The spin filtering stems from a combination 
of almost perfect matching of Gr and TM lattice con- 
stants and unique features of their electronic band struc- 
tures. Graphite films have occupied states at the Fermi 
level only around the K-point in the first (interface) BZ. 
Close-packed fee and hep Ni and Co have only minor- 
ity spin states in the vicinity of the same K point, at the 
Fermi energy. For a modest number of layers of graphite, 
transport from one TM electrode to the other can only 
occur via the graphite states close to the K point and 
perfect spin filtering occurs if the in-plane translational 
symmetry is preserved. For majority spins, the graphite 
film acts as a tunnel barrier while it is conducting for 
minority spin electrons, albeit with a small conductance. 

Compared to a conventional magnetic tunnel junction, 
a TM|Gr„|TM CPP junction has several important ad- 
vantages. Firstly, the lateral lattice mismatch is three 
times smaller than the 3.8% found for the now very popu- 
lar Fe|MgO|Fe(001) MTJs.^^ This will reduce the number 
of defects caused by strain that otherwise limits the thick- 
ness of the tunnel barrier and degrades the efficiency of 
spin injection. Secondly, the spin polarization approaches 
100% for an ideal junction with n > 3 graphene layers, 
and is only reduced to 70-90% for junctions with large in- 
terface roughness or disorder. Thirdly, the spin-filtering 
effect should not be very sensitive to temperature. From 
Fig. [SJa) and the corresponding figures for other thick- 
nesses of graphite, we see that the largest contribution 



to the majority spin conduction comes from tunneling 
at the M point where bulk Co and Ni have propagating 
states at the Fermi level and the distance in energy to 
states in graphite with the same k|| vector is a minimum. 
From Fig. [3] we see that the energy gap is almost 1 eV 
between the Fermi level and the closest graphite band at 
this point. To bridge the horizontal gap between states 
close to the K-point in graphite and the closest states 
in Co or Ni requires an in-plane momentum transfer of 
order Afc ~ 7r/a. The corresponding energy would be 
(comparable to) that of an optical phonon which is large 
because of the stiffness of a graphene sheet. 

To achieve perfect spin-injection into a single sheet of 
graphene is more troublesomc i^^'^^ The electronic struc- 
ture calculations presented here show that the carbon 
TT orbitals hybridize strongly with Ni (and also Co) sur- 
faces leading to the destruction of graphene's characteris- 
tic electronic structure. We have already suggested that 
dusting Ni (or Co) with Cu will lead to near-complete 
restoration of the graphene electronic structure because 
of the weak interaction between graphene and Cu. More- 
over Cu might also prevent rapid oxidation of the Ni(Co) 
(111) surfaces, which could be important for making 
practical devices. However, application of a bias would 
lead to a breaking of the translational symmetry respon- 
sible for the perfect spin filtering. The finite size of elec- 
trodes might also present a problem in practice especially 
if the potential drop occurs at the edges. The problem 
can be simply solved by forcing the electric field to be 
perpendicular to the TM| graphite interface as sketched 
in Fig. [TT] where the right electrode could equally well be 
placed on top of the graphite. 

In conclusion, we propose a new class of lattice- 
matched junctions, TM|Gr„|TM, that exhibit exception- 
ally high magnetoresistance effect which is robust with 
respect to interface disorder, roughness, and finite tem- 
peratures making them highly attractive for possible ap- 
plications in spintronic devices. 



11 



Acknowledgments 

This work is supported by "NanoNed", a nanotech- 
nology programme of the Dutch Ministry of Economic 
Affairs. It is part of the research programs of "Chemis- 
che Wetenschappen" (CW) and "Stichting voor Funda- 
menteel Onderzoek der Materie" (FOM) and the use of 



supercomputer facihtics was sponsored by the "Stichting 
Nationalc Computer Faciliteiten" (NCF), all financially 
supported by the "Nederlandse Organisatie voor Weten- 
schappelijk Onderzoek" (NWO). MZ wishes to acknowl- 
edge support from EU grant CARDEQ under contract 
IST-021285-2. 



^ V. M. Karpan, G. Giovannetti, P. A. Khomyakov, M. Ta- 
lanana, A. A. Starikov, M. Zwierzycki, J. van den Brink, 

G. Brocks, and P. J. Kelly, Phys. Rev. Lett. 99, 176602 
(2007). 

^ I. Zutic, J. Fabian, and S. D. Sarma, Rev. Mod. Phys. 76, 
323 (2004). 

3 A. H. C. Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, 
and A. K. Geim, Rev. Mod. Phys. -, (2008), cond- 
mat/0709.1163. 

* M. N. Baibich, J. M. Broto, A. Pert, F. Nguyen Van 
Dau, F. Pctroff, P. Etienne, G. Creuzet, A. Friederich, and 
J. Chazelas, Phys. Rev. Lett. 61, 2472 (1988). 

^ G. Binasch, P. Griinborg, F. Saurenbach, and W. Zinn, 

Phys. Rev. B 39, 4828 (1989). 
^ M. Julliere, Phys. Lett. A 54, 225 (1975). 
^ J. S. Moodera, L. R. Kinder, T. M. Wong, and R. Meser- 

vey Phys. Rev. Lett. 74, 3273 (1995). 

* T. Miyazaki and N. Tezuka, J. Magn. & Magn. Mater. 139, 
L231 (1995). 

® S. Yuasa, T. Nagahama, A. Fukushima, Y. Suzuki, and 

K. Ando, Nature Materials 3, 868 (2004). 
1° S. S. P. Parkin, C. Kaiser, A. Panchula, P. M. Rice, 
B. Hughes, M. Samant, and S. H. Yang, Nature Materi- 
als 3, 862 (2004). 

S. Yuasa, T. Katayama, T. Nagahama, A. Fukushima, 

H. Kubota, Y. Suzuki, and K. Ando, Appl. Phys. Lett. 
87, 222508 (2005). 

Y. M. Lee, J. Hayakawa, S. Ikeda, F. Matsukura, and 

H. Ohno, Appl. Phys. Lett. 90, 212507 (2007). 

S. Yuasa and D. D. Djayaprawira, J. Phys. D: Appl. Phys. 

40, R337 (2007). 
" P. X. Xu, V. M. Karpan, K. Xia, M. Zwierzycki, 

I. Marushchenko, and P. J. Kelly, Phys. Rev. B 73, 
180402 (R) (2006). 

^5 M. Zwierzycki, K. Xia, P. J. Kelly G. E. W. Bauer, and 

I. Turek, Phys. Rev. B 67, 092401 (2003). 
^'^ P. R. Wallace, Phys. Rev. 71, 622 (1947). 
" J. C. Slonczcwski and P. R. Weiss, Phys. Rev. 109, 272 

(1958). 

1® W. M. Lomer, Proc. R. Soc. London, Ser. A 227, 330 
(1955). 

1" T. Ando, J. Phys. Soc. Jpn. 74, 777 (2005). 

^° K. S. Novoselov, D. Jiang, F. Schedin, T. J. Booth, V. V. 

Khotkevich, S. V. Morozov, and A. K. Geim, Proc. Natl. 

Acad. Sci. U.S.A. 102, 10451 (2005). 

K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 

Y. Zhang, S. V. Dubonos, I. V. Grigorieva, and A. A. 

Firsov, Science 306, 666 (2004). 

K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 
M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, and 
A. A. Firsov, Nature 438, 197 (2005). 



2^ Y. B. Zhang, Y. W. Tan, H. L. Stormer, and P. Kim, 
Nature 438, 201 (2005). 

H. B. Heersche, P. Jarillo-Herrero, J. B. Oostinga, L. M. K. 

Vandersypen, and A. F. Morpurgo, Nature 446, 56 (2007). 

K. S. Novoselov, Z. Jiang, Y. Zhang, S. V. Morozov, H. L. 

Stormer, U. Zeitler, J. C. Maan, G. S. Boebinger, P. Kim, 

and A. K. Geim, Science 315, 1379 (2007). 
^'^ E. W. Hill, A. K. Geim, K. Novoselov, F. Schedin, and 

P. Blake, IEEE Trans. Mag. 42, 2694 (2006). 
■^'^ N. Tombros, C. Jozsa, M. Popinciuc, H. T. Jonkman, and 

B. J. van Wees, Nature 448, 571 (2007). 

H. Ibach and H. Liith, Solid-State Physics (Springer- 
Verlag, Berlin, Heidelberg, 1995), 2nd cd. 

Y. S. Dcdkov, M. Fonin, and C. Laubschat, Appl. Phys. 
Lett. 92, 052506 (2008). 
^° C. Oshima and A. Nagashima, J. Phys.: Condens. Matter. 
9, 1 (1997). 

Y. Gamo, A. Nagashima, M. Wakabayashi, M. Terai, and 

C. Oshima, Surface Science 374, 61 (1997). 

O. K. Andersen, O. Jcpscn, and D. Glotzel, in Highlights of 
Condensed Matter Theory, edited by F. Bassani, F. Fumi, 
and M. P. Tosi (North-Holland, Amsterdam, 1985), Inter- 
national School of Physics 'Enrico Fermi', Varenna, Italy,, 
pp. 59-176. 

P. E. Blochl, Phys. Rev. B 50, 17953 (1994). 
G. Krcssc and D. Joubcrt, Phys. Rev. B 59, 1758 (1999). 
G. Krcssc and J. Hafncr, Phys. Rev. B 47, 558 (1993). 
^® G. Kresse and J. FurthmuUer, Phys. Rev. B 54, 11169 
(1996). 

J. Neugebauer and M. Scheffler, Phys. Rev. B 46, 16067 

(1992). 

P. E. Blochl, O. Jepsen, and O. K. Andersen, Phys. Rev. 
B 49, 16223 (1994). 

G. Giovannetti, P. A. Khomyakov, G. Brocks, V. M. 

Karpan, J. van den Brink, and P. J. Kelly, Phys. Rev. 

Lett. 101, 026803 (2008). 
^° T. Ando, Phys. Rev. B 44, 8017 (1991). 

P. A. Khomyakov, G. Brocks, V. Karpan, M. Zwierzycki, 

and P. J. Kelly Phys. Rev. B 72, 035450 (2005). 
*2 K. Xia, P. J. Kelly G. E. W. Bauer, I. Turek, J. Ku- 

drnovsky, and V. Drchal, Phys. Rev. B 63, 064407 (2001). 
'^^ K. Xia, M. Zwierzycki, M. Talanana, P. J. Kelly, and 

G. E. W. Bauer, Phys. Rev. B 73, 064420 (2006). 

M. Zwierzycki, P. A. Khomyakov, A. A. Starikov, K. Xia, 

M. Talanana, P. X. Xu, V. M. Karpan, I. Marushchenko, 

I. Turek, G. E. W. Bauer, et al., phys. stat. sol. B 245, 
623 (2008). 

''^ M. Biittiker, Y. Imry, R. Landauer, and S. Pinhas, Phys. 

Rev. B 31, 6207 (1985). 
■'^ S. Datta, Electronic Transport in Mesoscopic Systems 

(Cambridge University Press, Cambridge, 1995). 



12 



I. Turek, V. DrchaL J. Kudrnovsky, M. Sob, and P. Wein- 
berger, Electronic Structure of Disordered Alloys, Surfaces 
and Interfaces (Kluwcr, Boston-London-Dordrecht, 1997). 
P. Soven, Phys. Rev. 156, 809 (1967). 
D. Glotzel, B. Segall, and O. K. Andersen, Sol. State 
Comm. 36, 403 (1980). 

T. Hahn, ed., International Tables for Crystallography, 
vol. A (Springer, 2005). 

G. Bertoni, L. Calmels, A. Altibelli, and V. Serin, Phys. 
Rev. B 71, 075402 (2005). 

To describe the adsorption of graphite on Co and Ni within 
model I we use two ESs with a WSR r = 1.12 a.u. at B 
and C sites, positioned 1.57 a.u. above the metal surface. 
For model II we use one ES with a WSR r = 1.7 a.u. at a C 
site, at 2.42 a.u. above the surface. If the distance between 
graphite and the metal surface is larger, which is the case 
for graphite on Cu, we use three ESs with WSR r = 1.6 
a.u., at A, B, and C sites, respectively, 3.33 a.u. above the 



surface. 

W. H. Butler, X.-G. Zhang, T. C. Schuhlicss, and J. M. 
MacLarcn, Phys. Rev. B 63, 054416 (2001). 
^■^ J. Mathon and A. Umerski, Phys. Rev. B 63, 220403(R) 
(2001). 

®* J.-C. Charlier, X. Gonze, and J.-P. Michenaud, Phys. Rev. 
B 43, 4579 (1991). 

G. Giovannotti, P. A. Khomyakov, G. Brocks, P. J. Kelly, 
and J. van den Brink, Phys. Rev. B 76, 073103 (2007). 
^'^ A. K. Geim and K. S. Novoselov, Nature Materials 6, 183 
(2007). 

®* D. J. Monsma, R. Vlutters, and J. C. Lodder, Science 281, 
407 (1998). 

^® K. Matsubara, K. Sugihara, and T. Tsuzuku, Phys. Rev. 

B 41, 969 (1990). 
®° G. Schmidt, D. Ferrand, L. W. Molenkamp, A. T. Filip, 

and B. J. van Wees, Phys. Rev. B 62, R4790 (2000). 



